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Recent studies on reproducing heat conduction in Hamilton systems revealed that anoma- 
lous thermal conductivity is ubiquitous in wide variety of systems even if those have strong 
mixing property, especially in those the total momentum is conserved. What is typical in 
those systems is that the thermal conductivity slowly diverges with system size while the tem- 
perature profile has well-scaled normal form. Those facts indicate that there are both diffusive 
and ballistic transport processes in such systems. Microscopic description and understanding 
for such non-diffusive energy transport, beyond its mere evaluation from the anomaly of the 
macroscopic transport coefficient, is expected so that it will be controlled and designed. The 
distribution of microscopic energy flux will be the most elementary quantity for this purpose. 
However, it is not clear what will be the normal distribution for microscopic energy flux since 
the energy flux has no direct relation with conserved quantities. In this note we discuss the 
distribution of microscopic energy flux in thermal equilibrium state. 

Let us consider the local energy flux J at a small volume element V in particle systems 
in equilibrium state. Microscopic expression of J is J = YliJi^ where j'j denotes the energy 
flux carried by each particle in the volume element V, hence the summation is taken for the 
particles in the volume. In the following we will focus on the ensemble distribution of j. 

For d— dimensional ideal gas, this local energy flux J is simply an average of single-particle 
energy fluxes and the each single-particle flux j is expressed as 

j^g.^^g.Z = (1) 

where e, v, p and m denote kinetic energy, velocity, momentum and mass of each particle, 
respectively. Since equilibrium distribution of momentum is 

P(p)dp = iVexp|^|dp {P = {kBT)-'), (2) 

where T, and N denote Boltzmann's constant, temperature, and a normalization factor, 



*E-mail address: shimada@ap.t.u-tokyo. ac.jp 



1/5 



J. Phys. Soc. Jpn. Short Note 

respectively, and the Jacobian for the transformation from p to j is calculated as 



dj 



2^) ^^\2^) ' 



dp 

the equilibrium distribution of energy flux is obtained as 



P{p)dp = Pip)^dj (4) 



Therefore 



J'0) = iV'j^exp|-(|)^/3ji|, (6) 

where N' is a normalization factor which does not depend on j. It is notable that the micro- 
scopic energy flux distribution in equilibrium state has a stretched exponential form, although 
it is just a transcription of Boltzmanian distribution function. 

For real gas with interaction between particles, the local energy flux J may not be deter- 
mined solely from the distribution of single particle flux since the correlation among them 
may be relevant. Yet the single-particle distribution is still essential and is a good observable 
for simulational physics. For the systems which have a Hamiltonian of 

i i<j 

the energy flux transported by a particle i can be expressed as 

where the first term is the same term as eq. (1). The second and third terms, come from the 
interaction, denote the advection term of the local potential energy and the energy transfer 
via the potential force, respectively. In the following, we will denote those three flux terms as 

3 = 3k + 3u + 3f with = SiPu3u = Ej and jp = - J2j {Qi " Qj) • • 

Since there is no correlation between the distributions of momenta and coordinates in 
the equilibrium ensemble, the expectation value of the flux for a given momentum p can be 
written as 

2 

< 3{P) >,= 1^ + P^jK + d, (9) 

where denotes the ensemble average over coordinate space, ^'(/3) is a tensor calculated 
as 

\ J I q 

and J =< jif + jp >q. Note that is essentially the virial. 
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Fig. 1. Distributions of microscopic energy flux. Points: P{j) obtained in Lenard- Jones particle sys- 
tem with e/3 = 1/3.4, and p ~ 0.25. Total number of particles is 843. Lines: theoretical predictions 
from the equation (6) and (11). 
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Fig. 2. Distributions of J/^, J[/, and Jp. Points: simulation results obtained from the same system 
as Figl. Lines: theoretical predictions. 



The equilibrium distribution of J is calculated as 

with a normahzation factor A^". This J, which stems from the potential and has p- first order 
form, becomes relevant for small p, and hence for small j. Since the flux j^j + j p is a product 
of p and potential-related term, the width of the distribution of itself is expected to be broader 
than the one of J and they are related by 

^Og{(yju+jF ) = log(crvi,(^) ) + log(crj) . (12) 

Eqs.(6), (9), and (11) imply that the distribution of the total local energy flux P{j) shows 
crossover from Arrhenius-type region, in which potential term is dominant, to non-Arrhenius- 
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type region, in which momentum energy advection is dominant. 

We now examine the argument using three dimensional Lenard-Jones particle system with 
the potential 



and the uniform mass (m = 1). This system is known to have normal thermal conductivity.^ 
The distribution P{j) is estimated by computer simulation at high temperature (e/3 = 1/3.4) 
and not dense (number density p = 0.25) condition (Fig.l). We can sec that the distribution 
with crossover discussed in this note appears in the simulation result. Since the density is low 
and the temperature is high, it is expected that ^{(3) can be well estimated by taking the 
average for single pair using two-body distribution. The lines in Fig.l show the distributions 
calculated using the analytic result for the second virial coefficient.^ As shown in Fig.l, the 
envelope of theoretical predictions for P(jx), P{ju), and P{jF) well explains the profile of 
simulational result of P{j) i.e. the peak position and the stretched exponential tail. It is also 
confirmed that theoretical prediction and simulational result show good agreement in each 
P(ji^), P{ju), and P{jF), except for large jp regime where the latter has stretched exponential 
tail which originates from high-energy collision events (Fig. 2). 

We have focused directly on the distribution of the microscopic energy flux. From the 
microscopic point of view, there arc many ways to transport the energy without average 
momentum flow. Asymmetry in the higher order cumulant of P{p) is a possible candidate, 
although it is known to be usually hard to detect the deviation directly in the distribution 
of dynamical variable. Correlations between momentum variables and coordinate variables 
will be an another possibility. For local heat flux, spatial correlations among particles is also 
essential for anomalous behaviors and that should be detected as the deviation from the 
independent sum of the single particle distribution. Furthermore, to consider the boundary 
in non-equilibrium states, for instance,'' anomaly can be estimated only from microscopic 
quantities. The argument we took and the obtained normal distribution give good basis for 
considering those problems. 
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